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ABSTRACT 

This paper provides two families of flexible and simple galaxy models. Many representatives of these fami- 
lies possess important cosmological cusps, with the density behaving like r _1 , r~ 4 ^ 3 , or r~ 3 ^ 2 at small radii. The 
density falls off between r~ 3 and r~ 5 at large radii. We provide analytic and anisotropic distribution functions 
for all the models. Unlike many existing methods, our algorithm can yield tangentially anisotropic velocity 
dispersions in the outer parts, and so is useful for modeling populations of satellite galaxies and substructure 
in host galaxy halos. As an application, we demonstrate the degeneracy between mass and anisotropy for the 
satellite galaxy population of the Milky Way. This can introduce a factor of ~3 uncertainty in the mass of the 
Milky Way as inferred from the kinematics of the satellite population. 

Subject headings: galaxies: kinematics and dynamics- methods: analytical - stellar dynamics 



1. INTRODUCTION 

One of Eddington's famous discoveries is that the is otropic distribution fun ction (D F) of a sphe rical stellar system can be 
calculated from the density using an Abel transform pair (Eddingto nTl916tlBinnev & Tremainel 19871 p. 237): 



f(E) 
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(1) 



Here, E — if/ - v 2 /2 is the binding energy per unit mass, and if/ is the relative potential. 

However, galaxy hal os produced in cosmological simulations are held up by anisotropic velocity disper sions (see e.g., 
Hansen & Moore 2006). The recovery of an anisotropic DF for a spherical system is a more difficult problem ( Hunter & Oian 
1993). Now, the steady state DF can depend not only on E but also on the magnitude of the specific angular momentum L — rvj 
through Jeans' theorem. Here vt is the two-dimensional velocity component projected on the tangential plane. 

At least in the comp arison of observational data with models (e.g., Ivan d er Mareletal T l2000t iWilkinson et ail l200l 
iTreu & Koopmansl2004l) . a widely used Ansatz for the form of DF is 

f{E,L) = L- 2 Pf E {E), (2) 

where /e(E) is a function of the binding energy alone. It is easy to show that the models generated by this DF exhibit the property 
that 



0=1- 



(v 2 T ) 
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(3) 



Here, (v 2 ) and ( y 2 ,) are the r adial and the tan gentia l velocity second mo ments. The general inversion formula for the unknown 
function f E (E) jCuddefordlU99ll iKochanekll 19961 IWilkinson & Evansll 19991 see also Appendix |Aj is no more difficult than 
Eddington's formula (eq.[Q. 

We note that it is common practice to characterize the velocity anisotropy of a spherically sym metric stellar system by th e 
parameter defined in equation 0, which is sometimes referred to as the anisotropy parameter (Binnev & Tremaine 1987). 
If < /3 < 1, the velocity dispersion ellipsoid is a prolate spheroid with the major axis along the radial direction (radially 
anisotropic), while it is an oblate spheroid with the tangential plane being the plane of symmetry (tangentially anisotropic) if 
P < 0. In particular, ft = 1 and (3 = -oo indicate that every star is in a radial or circular orbit, respectively. On the other hand, 
systems with isotropic velocity dispersions have /3 = 0. In general, the anisotropy parameter varies radially, but for the system 
with a DF of the form of equation 0, it is constant everywhere. 

Despite the attractiveness of a simple inversion of the DF in equation ©, the assumption of uniform anisotropy is not always 
desirable. For example, it is useful to be able to build mod els that are isotrop ic in the center and radially anisotropic in the outer 
parts, as seems to be the case for dark matter halos (Hansen & Moore 2006). DFs that are tangentially anisotropic in the outer 
parts are helpful for understanding the substructure and the satellite galaxy populations in dark matter halos. 

In this paper we study algorithms for building galaxies with varying anisotropy, using sums of DFs of the form of equation©. 
§ |2] introduces a new family of cusped halo models, the generalized isochrones, and § [3] examines the generalized Plummer 
models. Every member of the family has analytic potential-density pairs, velocity dispersions, and DFs. Specific examples 
are provided with the cosmologically important r~ l , r~ 3 ^ 2 , and r~ 4 ^ density cusps. §0]discusses a specific application of our 
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Fig. 1. — Density profile of the generalized isochrone models. The classical isochrone (p = 2) is cored, whereas all other models are cusped with p ~ r * p 2 * at 
small radii. The small arrows indicate the half-mass radius. 



theoretical models to the problem of estimating the mass of a host galaxy from the motion of its satellites. Finally, §[5]concludes 
with a brief comparison between our methods of producing varying-anisotropy DFs and other suggestions in the literature (e.g., 
lOsinkovll 979tlMerrittn 985tlCuddefordll 99lHGerh a rdl1 9911). 

2. THE GENERALIZED ISOCHRONE SPHERES 
The isochrone sphere has the potential-density pair (Henon 1959; Eg gen. Lvnden-Bell. & Sandaii IT962t 



Evan s, de Zeeuw. & Lvnden-Belll l 990 ) 



GM 
b + s ' 



P = 



M (2s + b)b 
4n s 3 (b + s) 2 ' 



(4) 



where M is the total mass of the system, s - Vfr 2 + r 2 , and b is a constant defining the scale length. Like the Plummet] ( 1191 It) 
model, the density is finite and well behaved at the center. Unlike the Plummer model, the density falls off as r at large radii, 
which is a better fit to elliptic al galaxies than t he r fall-off of the Plummer model. There are many available studies regarding 
its DFs (e.g., Gerha rdTl991tlBertin et all 19971) . some of which are essentially numerical in nature. 
Inspired by this, let us consider the potential-density pair 



GM 

b + s' 



M2brP + bP(\ + p)(b + s) 
\n r 2 -Ps 2 P- l (b + s) 3 



(5) 



where 



s = (b p + r p ) 



p\ l /p 



(6) 

Here, p > is a new parameter used for the generalization of the isochrone sphere (eq. 0}. The condition that the density 
decreases monotonically with increasing radius restricts p < 2, with the limit (p = 2) being the cored classical isochrone sphere. 
In fact, all the halo models with < p < 2 are cusped. In particular, the density behaves as p ~ r~ {2 ~ p) as r — > 0. On the other 
hand, we find that p ~ r~ q where q = min(4, p + 3) as r — > oo. The density profiles of some members of this family are shown in 
Figure [2 We note that the case p = 1, for which ift = GM/{2b + r), is actually recognized as the Hernquist ( 1990) model, albeit 
with a slightly unusual identifica tion fo r th e scale length (a = 2b). Models with p — 1/2 and p = 2/3 possess the cosmological 
cusps suggested by Moor e et ail (119981) and lEvans & CollettJ ill 9971) . 
The enclosed mass and the circular speed are found to be 
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(7) 



The potential ip(r) of the generalized isochrone sphere is simple enough that the inverse, r(i[r), can be easily found (with G 
b = 1 without any loss of generality): 



M : 



m = 



(8) 



Using this, in principle, one can construct infinitel y many di fferent expressions for p(i[r, r), each of which can be inverted to yield 
the DF by means of integral transforms (Lvnden-Bell 1962). However, for most cases this procedure is analytically intractable. 
Rather, it is often more productive to look for a specific form of p(i//, r) that is associated with a plausible Ansatz for the DF, 
whose inversion reduces to elementary integrals. 
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2.1. Two Component Distribution Functions 
One such possible (r,t/r)-split with an analytically tractable inversion is given by 

p = Ell r"- 2 if, 2p+l (l - if,) 1 - Zp + — r 2 "-V 2p+2 (l - <A)'~ 2p - (9) 
4n 2n 

We note that a DF of the form of equation (0 can be inverted from an (r,i/')-split that is a monomial of r, with its power index 
linearly related to ft. By having the (r,(/»)-split given by a linear combination of two monomials of r, the inversion of equation 
can be achieved by simply extending the known procedure, while the system is no longer restricted to have uniform velocity 
anisotropy. Inspired by this, let us suppose that the DF is of the form 

f(E, L) = LP' 2 ME) + L 2 P- 2 f 2 (E), (10) 

where f\(E) and f%(E) are functions of E to be determined. Then, the density is found to be 

P F0V2+1/2) Jo W Jl T(p +1/2) Jo W 

where F(x) is the gamma function. By comparing this to equation a possible DF may be found by inverting the integral 
equations 

Jo 2^ 2+5 /V/ 2 r(/?/2) r ^ 

JT*E (* - Ey~ 1/z f 2 (E) = ^ +2 d - (12) 

The inversion procedure is essentially identical to the one that leads to equations iA2\ and jA3> . After this inversion, we find that 
fi(E) and /2(F) are 

- ^ + *, + (13, 

both of which are always nonnegative for all accessible ^-values (0 < E < if/ < 1 /2). While we have given the result in complete 
generality in terms of hypergeometric functions 2^1(0, b; c; x), the DFs reduce to entirely elementary and analytic functions if 
p = 1/2, 1, or 2. In fact, if p = 1/2, both hypergeometric functions are just unity, whereas if p is an integer, both reduce to 
elementary functions. If p is a half integer other than 1 /2, only f 2 reduces to an elementary function. 

Of particular interest is the case in which p = 1 /2, for which the model possesses a central cusp like p ~ r~ 3 ^ 2 , as found by 
cosmological simulations (Mo ore et al.ll998 p. The corresponding DF is very simple, namely, 

3 E 2 2 5 ' 4 ■ 3 E 5 ' 4 
/(£ ' L) -^T + 57r 5 /2F(l/4) 2 Z^- (14) 



If p = 1 (the Hernquist model with a rescaled scale length), 



f(E,L)= ' 



2 7 /V(l -E) 2 
or if p = 2 (the classical isochrone sphere); 
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(16) 



2 15 / 2 tt 3 (1 -Ef 

Both of these DFs are no more complicated than the isotropic DFs deduced by Hernauist ( 1990) and H enonl i ll 9601) . respectively. 

2.1.1. Kinematics 

To find the velocity dispersion of these models, we exploit the splitting p — p\ + p 2 , where 

f „o , (p + l)rP- 2 il/ 2 P +l (l-tl/) l - 2 P C ,„ , , r 2 P- 2 i]/ 2 P +2 {\-ib) l - 2 P 

Pt=J LP- 2 ME)d\ = — ^± W - ; p 2 = J L 2 P- 2 f 2 {E)d\ = v ^ V) . (17) 

Then, 

p<v 2 > = ^-y, + J^-V 2 ; P <v 2 ) = J^V, + ^ V 2 , (18) 

2(/? + 1) 2p + 3 2{p + 1) 2p + 3 
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Fig. 2. — Behavior of the anisotropy parameter generated by the DFs (eg. 1101 for the generalized isochrone spheres. From top to bottom, the parameter p varies 
from 1/4 to 2 with increments of 1/4. Note that the last curve, p = 2, corresponds to the classical cored isochrone sphere. 



where 



Vi = (1 -#) 2 Fi(3 + i, \;2p + 2 + i;iff) = zFA2p- 1, l;2p + 2 + /; 



The anisotropy parameter 



2 



(2p + 3)Vi + 8rPifrV 2 



(19) 



(20) 



(2/? + 3) Vi +4rPtl/V 2 

is now no longer constant. Figure [3] shows the behavior of ft for some of these models. The curves can be understood on 
examination of equation J20I . This reveals that whereas /? = 1 - (p/2) at r — for any p, lim^oo/? = 1 - p if p > 1, or 
lim r ^oo/? = 1 - (p/2) if p < 1. At the critical value p = 1, we find that the anisotropy parameter monotonically decreases from 
[i — 1/2 at r = to p 1 — >5/18asr— > oo. Hence, the DFs discussed in this section build the generalized isochrone spheres with 
either decreasing (3 in the inner region and increasing ft in the outer region (p < 1), or decreasing p 1 throughout (p > 1). The 
particular case of the latter also includes the Hernquist model (p = 1) and the classical isochrone sphere (p = 2). 
The simplest case is again p = 1/2, which leads to 

2 + 3r 1/2 + r . 2p\ + 3p?. \jj 2 + 4r l/2 + r 



<K) = 



4pi + 3p2 4> 



■V, 



pi+pi 12 6 + Wr 1 ' 2 + 3r 
Then we find that the virial theorem holds locally for this case: 



<v4> = 



(vj) + <V 2 T ) = 



and the anisotropy parameter 



_ 1 6 + 8r 1/2 + 3r 
^"4 2 + 3r'/ 2 + r 
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6pi + 6p2 4 1 4 1 
p\+ pi 12 2 ' 



3(r-2) 2 -r 1/2 (r + 2) 



(21) 



(22) 



(23) 



4(r-4)(r- 1) 

varies from ^0 = 3/4 at r — 0, tojS = 1/ V2 at r - 2, and back to /3 — > 3/4 as r — > oo. This model in fact corresponds to one of the 
generalized hypervirial models studied bv lAn & EvansH200 5a). That is, ij/ — [1 + (1 + r 1 ' 2 ) 2 ] -1 = (2 + 2r'^ 2 + r) _1 so that ro = 2fe 
and (p, c) — ( 1 , 1 / V2) in their parameters. 

We note that both V\ and V 2 in fact reduce to the incomplete Beta functions, which may be useful for the purpose of numerical 
calculations. Furthermore, if 4p is an integer, both Vi and V 2 reduce to elementary functions of ifr, and so do the velocity second 
moments. In particular, if p — 1 [for which iff = (2 + 

128 + 303r + 250r 2 + 90r 3 + 12r 4 

; (24) 



(vt) = ~(l+r)(2 + ry 
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<v 2 >= -(1 +2r)(2 + r) 3 In 
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24(2 + r) 

64 + 217r + 21 lr 2 + 84r 3 + 12r 4 
12(2 + r) ' 



(25) 
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3. GENERALIZED PLUMMER MODELS 
Next, let us consider the one-parameter family of potential-density pairs, 

GM = (p + 1)M gP 

(aP + rPy/p' 4-n r 2 -P(aP + rP) 2+1 /p' 

which was originally introduced by Veltman nl ill 9791). This includes the Hernquist ( 1990) model (p — 1) and the Plummet] d!9 1 ll) 
model (p = 2) as particular cases. Evans & An (2005) recently found that this family can be constructed from power-law DFs of 
the form of f(E, L) cc IP- 2 E® P+V> I 2 . Here we construct slightly more complicated DFs that build this family of models. We first 
find every DF of the form of equation (0, and then use them to build models with a more general variation of anisotropy. 

3.1. Distribution Functions with Constant Anisotropy Parameter 
With G = M = a = 1, we find that 

p-Pp = lll^-^i _ ^i-VX-Wp (27) 
4tt 



and that 



for m < p + 3 — 2/5. Then, equation iA2\ reduces to 



= (28) 

i/<=o 



HE n = 2 ^ P + 1 1 d " +1 f E d A ( 79) 

n ' ' (2tt) 5 /2 F(« -j8- 1/2)F(1 -P) L 2 P dE" +l J (1 - ipP) A (E - ^V 2 ^-" ' ^ ; 

where « = L3/2 — /3J is the integer floor of 3/2 - B, and A + 1 = 2( 1 - p)jp. If A > 0, it is at least formally possible to derive the 
series expression of the DF for arbitrary p and j3 from the integral form, namely, 

f( ' (2tt) 5 / 2 T(l -y8) L 2 ! 3 f-J T{pk + p + 5/2-/3) T{A) k\ ' 

which becomes the generalized hypergeometric series if p is a rational number. On the other hand, if -1 < A < 0, we find 
that lim^i- f(E,L) < so that the corresponding DF is unphysical. Hence, for the potential-density pair of equ ation (1261 . the 
constant anisotropy DF is physical only if p < 2(1 - /3). This means that the hypervirial models of IE vans & Anl (12005). which 
have an anisotropy parameter (3 — 1 - (p/2), have the maximally radially biased velocity dispersions for a given p and a constant 
B. Although there exist models with /3 locally exceeding 1 - (p/2) for the generalized Plummer sphere of given p (see e.g., 
Baes & Deionghe 2002), /? at the center cannot b e greater than 1 - (p/ 2) for physically valid DFs. This is in agreement with the 
cusp slope-central anisotropy theorem derived bv lAn & Evans! (2005b). 

The corresponding velocity dispersions can be found either from equation (IA8I or by solving Jeans equation. The results can be 
generally expressed using the incomplete Beta functions. The equivalent results expressed in terms of hypergeometric functions 
are also found in equation (9) of Evans & An] ( 120051) . 

For particular values of p or /?, the DF of equation ( I29l l (or equivalently, eq. I30l > reduces to a simpler form. For example, if 
(3 = 1 /2, the DF is (see eq. |A3> 

p + 1 EP +1 
(2nfL (1 - EP) 1 Ip 



f(E,L) = ^ - ^, ;n [{p + 2) - (2p + l)EP] , (31) 



where the positive definiteness of the DF implies that < p < 1. Similarly, for p = -1/2, the DF is 

f(E,L) = (27r) 3 JiZe-pJWp i (p + 3){p + 4) " (5/ ' 2 + l2p + 3W + 2p(2p + 1)e2 "\ ' (32) 

which is nonnegative for < E < 1 if < p < 3. On the othe r hand, the DF of the Hernquist model with constant [} reduces to 
the hypergeometric function (see e.g. JBaes & Deionghel2002l) 

2^F(5-2/3) E 512 -! 3 I 7 \ 

f(E,L) = tt, 1 — ^-2^1 1 -28,5 - 28;- -8,E\, (33) 

J ' (2nfl 2 Y(l 12 - B)Y(\ - (5) L 2 P \ H P 2 P J ' 

which is nonnegative everywhere if /3 < 1 /2, whereas that for the constant-/? Plummer model reduces to 

3 • 2^'F(6 - 2/3) E 112 - 
(2w) 5 / 2 F(9/2-/3)r(l -p) ~L 2 P 



3 • 2^'r(6 - 2/3) E 112 -? ( 1 9 - 28 11 -28 , , 



which is physical if B < 0. 
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Fig. 3. — Behavior of the anisotropy parameter generated by the DFs (ea. 1421 for the generalized Plummer models. Solid lines, p = 1 (Hernquist model); dotted 
lines, p = 3/4; short-dashed lines, p = 1/2; long-dashed lines, p = 1/4. 

3.2. Distribution Functions with Outwardly Decreasing Anisotropy Parameter 

Since equation d27> is valid for any constant p with common p, it is, in general, also possible to express the density profile of 
the generalized Plummer model by 



-2(l-ft)/p . 



4tt 



2> - 1. 



(35) 



where the parameters p, > are normalized weights. Consequently, the superposition of the constant anisotropy DFs, f(E, L) = 
2, Pifi(E, L), still builds the generalized Plummer sphere. Here each fi(E, L) is given by equation ( I29> (or equivalently, eq.!30> 
in which p is replaced by /?,-. The velocity second moments for these models are found to be 

2 v tiiS(r;p,Pi) 2 v 2p.j(l -Pi)S(r;p,Pj) 

P (v r ) = P if, ^ p + 4 _ 2/? . ; p<V = p-A 2j — ~ -n — ■ (36) 



p + A-lpt 



where 



SU, /,,) = A 1 -^- ~ U; ^ + 2; -1) = W + 3, 1; 1^ + 2; , 



(37) 



p p rP 

and therefore the anisotropy parameter /? is no longer constant. 

For the particular case in which the sum only contains two terms (such that p\ — p., p2 — 1 an d Pi - for all other i values) 
the anisotropy parameter is 

p lM (p + 4 - 2p 2 ) S(r; p,fi{) + p 2 (l - p)(p +4-2/30 S(r; p,p 2 ) 



P = 



(38) 



p( P + 4 - 2p 2 ) S(n P ,Pi) + (l - p)( P + 4 - ipo S(n P ,p 2 ) 

which is a decreasing function of r, provided that < p < 1 and P\ + p 2 . If we assume P2 < Pi < 1 - (p/2) and < p. < 1, the 
limiting values are found to be 

Ptpip + 4 - 2,62) +J S 2 (1 - p)(p + 4 - 2ySi) 



limyg = 



On the other hand, 



p{p + 4 - 2j&) + (1 - R>(p +4-2A) 
P(r = 0)=pi if 



^ r = 0) = n own ^1 ^l - if /?i < 1 - p. 

The simplest example of the DFs of this kind is obtained when we choose P\ = 1/2 and /% = -1/2: 



(39) 
(40) 
(41) 



f(E,L) 



(2^)3 



(p + 2)-(2p + l)£P | 2 (g + 3)Q> + 4) - (5p 2 + 12/7 + 3)E" + 2p(2/> + 1)£ 2 ^ 1 

which is physical if p < 1 . The behavior of the anisotropy parameter p for the model given by DFs of equation (1421 is shown 
in Figure |3 Again we note that equation d37t in general reduces to the incomplete Beta function while, for this case, it further 
reduces to an expression involving only elementary functions of i[> if 2p is an integer, and so and so do p and the velocity second 
moments. 
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Fig. 4. — Curves illustrating the degeneracy between mass and anisotropy. All the models have the same half-mass radius and reproduce the same mean kinetic 
energy (T rr ) of the satellite population of the Milky Way. The inferred mass of the Milky Way (in units of 10 12 M Q ) is plotted against the parameter p, which 
controls the velocity anisotropy. 



As an application of the preceding theory, let us consider the problem of estimating the mass of the Milky Way from the radial 
velocities of its satellites. There are 27 distant globular clusters and dwarf galaxies at Galactocentric distance r greater than 20 
kpc (see e.g., tables 2 and 3 of Wilkinson & Evans 1999). These range from Arp 2 at r = 20 kpc out to Leo I at r = 254 kpc. The 
median distance of the satellite population is r me d « 38 kpc. The num ber densi ty of the s atellites falls off roughly like r -3, . This 
dataset is often used to estimate the mass of the Milky Way (Kochanek 1996; Wilkinson & Evans 1999), as the H I rotation curve 
cannot be traced beyond ~20 kpc. The satellites are therefore the only available probe of the gravitational field of the distant halo. 

All the satellites have accurate radial velocities, but proper motions are available only for a handful of the closest objects. 
Of course, the heliocentric radial velocity and the Galactocentric radial velocity are basically the same (once the motion of the 
Sun has been taken out), as the Galactocentric distance of the Sun is much smaller than those of the satellites. The kinematic 
data therefore constrains the mean kinetic energy in the radial direction T rr = (v 2 ). Again using the data in tables 2 and 3 of 
Wil kinson & Evans! d!999l) . we find that T rr * ( 1 29 km s" 1 ) 2 . 

To model the Milky Way halo, we assume that its potential and density take the form of either the generalized isochrone 
or th e gene ralized Pl ummer models. Furthermo re, we assume that the number density of satellites shadows the total density 
(e.g., Kochanek 1996; Wilkinson & Evans 1999). This means that the satellites are drawn from DFs of the form given by either 
equation dlOt or equation J42l >. We set the scale length so that the half-mass radius of the assumed density profile corresponds to 
the observed median radius of the population. The mass of the Milky Way M enters into the potential and the DF. It can therefore 
be constrained by requiring that the mean kinetic energy in the radial direction be consistent with the data, viz, 



The parameter p in the DFs controls the variation of the velocity anisotropy with radius. Hence, equation J43i defines a line in 
the plane (M, p), which is the mass-anisotropy degeneracy curve for the model. Figure |4] inferred from numerical integrations, 
shows such curves computed for a number of the generalized isochrones and the generalized Plummer models. 

In this problem, the mean kinetic energy in the radial direction is fixed. Models with tangential anisotropy at large radii therefore 
lead to a higher inferred mass for the Milky Way galaxy. The relative location of the curves in Figure |4] can be understood by 
examining the behavior of the anisotropy parameter for the models. More importantly, as the velocity anisotropy changes, the 
inferred mass of the Milky Way ranges from 0.8 x 10 12 M Q to 2.5 x 10 12 M Q . Now, the velocity anisotropy of the satellite galaxy 
population is unknown, so the mass-anisotropy degeneracy by itself enforces a factor of ~3 uncertainty in the mass of the Milky 
Way. This fundamental limitation can only be overcome by proper-motion data. 

Of course, there are other uncertainties; for example, the half-mass radius of the satellite population used to set the scale 
length in these calculations is also not well known (note that the locations of the satellites are likely to be biased toward their 
apogalacticon) and has an uncertainty of ~20%. This carries over into an additional ~20% uncertainty in the inferred mass. This 
is small in comparison with the mass-anisotropy degeneracy, but not insignificantly so. 



4. AN APPLICATION: THE MASS OF THE MILKY WAY 



Jp{v 2 r )d 3 r 



* 16,640 km 2 s" 2 . 



(43) 



5. CONCLUSIONS 



This paper provides a number of flexible and simple galaxy models extended from two classical cored profiles: the isochrone 
sphere and Plummer model. We also derive analytic distribution functions (DFs) with varying velocity anisotropy for all the 
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models including new DFs for the classical isochrone sphere and the Hernquist model. 

DFs of the form f(E, L) = L^/e{E) build a spherical system with a constant anisotropy parameter ft, where E is the binding 
energy and L is the specific angular momentum. The assumption of constant anisotropy is not enough to provide realistic models, 
as both observational data and numerical simulation s suggest that dark halos a nd their co nstituents (such as satellites) have 
velocity anisotropics that vary with radius (e.g., White 1985; van derMareI199HlHansen & MoorefcOOfl) . 

Here we have explored models whose DFs are superpositions of two or more such terms, namely, 

f(E,L)=Y i L- 2 *ME). (44) 

i 

This can be done in two ways. First, as illustrated by the generalized isochrones in §[2] we can sometimes find a splitting of the 
density in the form 

P(r) = J] r- 2Pi gi(ifr(r)). (45) 

i 

and apply the inversion separately to each component. For the generalized isochrones, this provides models with either decreasing 
anisotropy parameter ft in the inner region and increasing ft in the outer parts, or decreasing ft throughout. Second, as illustrated 
by the generalized Plummer models in § [3] we can always add together weighted sums of constant anisotropy DFs, each of 
which individually reproduces the required density. This appears to provide models with decreasing anisotropy parameter ft 
throughout. At large radii, the DF is dominated by the component of smallest ft, and the models can be either radially or 
tangentially anisotropic in the outer parts, according to choice. For both cases, it appears that the anisotropy parameter always 
decreases on beginning to move outward from the very center. 

There have been a number of previous algorithms suggested for building anisotropic DFs with varying anisotropy. For example, 
the Osipkov-Merritt al gorithm ( Osi pkovl 1 979tlMerrittl 1985) builds models that are isotropic in the inner parts and tend to extreme 
radial anisotropy (ft — > 1) in the outer parts. Gerhard (1991) provided a clever algorithm for building separable DFs of the form 

/ = g(E)h(L/L c (E)), (46) 

where L c is the angular momentum of a circular orbit with energy E. The analytic examples provided by Gerhard i 19911) all have 
increasing radial anisotropy in the outer parts. His method can be used to generate tangentially anisotropic models, but usually at 
the cost of a numerical inversion. The methods developed in this paper therefore provide a complement to the existing inversions. 

Tangential anisotropy often arises in accreted populations, such as satellite galaxies. As an application of our models, we have 
shown that the mass-anisotropy degeneracy by itself provides a factor ~3 uncertainty in the mass of the Milky Way galaxy as 
deduced from the kinematics of its satellites. 

The density profiles of models studied in this paper typically fall off faster than p ~ r~ 3 and slower than p - r~ 5 at large radii, 
which are in good agreement with both observations of spheroidal components of galaxies and results from numerical simulations. 
At small radii, they are typically cusped, with cusp indices in the range suggested by numerical simulations. In particular, there 
are mem bers with density cusps like p ~ r" 1 , ~ r~ 4 ^ 3 , and ~ r~ 3 ^ 2 , which have b een suggested as important on cosmogonic 
grounds (Navarro. Frenk. & WMtd[l995tlEva ns & Co llettll l997: Moore et al. 1998), althou gh observational evidence seems to 
favor the presence of a c onstant-densit y core (iTvson. Kochanski. & deirAntonioll998tlPalunas & Wilhamsll2000Ude Blok et alJ 
1200 It iKlevna et alJl2003l iDonato. Gentile. & Saluccill2004) . However, even if real galaxies are cored, the core would be such a 
small fraction of the halo that it is not practically important for studies of the whole. Our families of new models therefore should 
find widespread application in the modeling of galaxies and dark halos. 

We thank C. Hunter, who made a number of interesting comments on the various incarnations of this paper. 
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APPENDIX 

A. CONSTANT ANISOTROPY DISTRIBUTION FUNCTIONS 

In this appendix we summarize some results regarding constant-anisotropy DFs used in the main body of the paper. For the DF 
given in equation Q, direct integration over velocity space gives 



H 7PTOI2-P) J,, W ' J y ' ' 

The formula for inversion for the unknown function /e(E) is given by (Cuddeford 1991 ) 

1 2? d C E dtfi d"h 



1}P (2n?l 2 T(\ -a)T{\ - p) 



— - f 

~P) dE J,, 

I 



f(E,L): 

w j yLTiy-v (L — a)i {L — p) an, j 1° ~ WF a W" 

1 TP r C E dtff d n+l h 1 d n h 

+ — — , (A2) 



(E - tff) a dif/ n+x E a diff" 



!//=() 



where h{ifr) = r 2/3 p is expressed as a function of \ji, and n = [3/2 — p\ and a — (3/2 — /?) - n are the integer floor and the fractional 
part of 3/2 — p, respectively. This includes equation as a particular case of ft — (n = 1 & a — 1/2). If p 1 is a half-integer 
constant (i.e.,/? = 1/2, -1/2, -3/2, and so on), the above inversion simplifies to 



1 1 d 3 ' 2 ^h 

f(E,L) = 



2n 2 L 2 P(-2p)\\ d^l 2 ~P 



(A3) 

ip=E 



that is, this only involves differentiations in the process ( Cudd eford! 1 99 11) . In addition, the expression for the differential energy 
distribution (DED) is found to be 



dp 
dE 



. - J -jp-W-Eo)* = m 3/2-f» ^ ME )m-E y, (A4) 

dM r dp , (2^) 5/2 F(l -P) C E mo r on m 



where >J/(rE) = E, and 5{x) and ®(x) are the Dirac delta 'function' and the Heaviside unit step function, respectively. In particular, 
if p = 1/2, the last integral does not explicitly involve the potential, and so the DED may be found simply as dM/dE = 
4n 3 r 2 f E (E). 

It is also straightforward to find the integral for the second velocity moments; 



p{v 2 } = (2 f ' 2 F(1 ® \ ty-E) 3 > 2 -ef E (E)dE; (A6) 



0A- 

o 



P<Vt>= "^ r ' (5 : 2 " m I (ifr-EY'^ME)dE = 2(1 -fifty). (A7) 

In fact, the velocity d ispersio n of any co nstant anisotropy models may be f ound from h(tf/) without explicitly evaluating the DF 
dDeionghell987UBaes & De ionghe 2002). On differentiating equation ( IA6> with respect to ift and noting that the right-hand side 
is the same as equation JA1I . that is, d{p{v 2 )) I d\f/ = p, it follows that 

i frf i r>l> \ C d\li' 

?> = - di// p(r,f) = — dfhm = ^ r dr'Jjr'^p. (A8) 
P Jo «W Jo ^P Joo dr' 

We note that equation (IA8I can also be derived fro m the integral solution of Jeans' equation with a constant p (see eq. 9 of 
IE vans & Anll2005l) . This also indicates that equation dA8l > is actually a general result for a constant p model, regardless of the 
specific form of the DF that generates the model. 

As an example, let us think of the DFs that build the generalized isochrone sphere with a constant anis otropy parameter. First, 
the DF and the DED for the generalized isochrone sphere with p — 1/2 can be found from equation ( IA3t : 

1 E 2 g(E) dM [(1 - E)P - E' } ] 1/p 

f(E L) — ' = — - 2(E) (A9) 

J 7 ' (27r)iL(l-E) 2 n[(l-E)P-EP] llp ' dE 2(1 - E) 2 p SK 

g(E) = E 2p - 1 [6(p + l)E - 6E 2 - (p + l)(2p + 1)] + (1 - E) P E P ~ X [(p + l)(p + 2) + 12£ 2 - 6(p + 2)e\ + 6(1 - E) 2p+l . (A10) 

Here, because < E < ifr < 1/2, the nonnegativity of the DF is satisfied only if < p < 1. In addition, the integration involved 
in the inversion (eq. lA2> is also analytically tractable if p = 1 - (p/2). That is, for the density profile of the generalized isochrone 
sphere, we find that 

r 2 - p p = ^" +2 d ~ *) l - p + ~ ^ 2P+2(1 ~ ' (AU) 



2(2nf 12 T{2-P) 



PL r 

-p)Jo 
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This leads to the DF of the form of equation (0 
1 1 



f(E,L)-- 



[2-p 2P/2+3/2 7r 5/2 r( - ; , / 2) 



(A12) 



which build a constant anisotropy model with /3 = 1 - (p/2). 
If /? = 1 (corresponding to the Hernquist model) and /? = 1/2, the DF and the DED further reduce to a particularly simple form; 



f(E,L) = 



3 E 2 



dM 



6(1 - 2E f 



(A13) 



2?r 3 L ' dE 

which was originally found bv lBaes &Deionghdll2()0l for the lHernauistll l990> mode l. For p — 2 an d /3 — 0, the well-known 
expression for the DF of the isotropic isochrone sphere (see e.g., Henon 1960; Binnev & Tremaine 1987, p. 239) is recovered; 

1 



/(£) = 



215/2^3(1 _£)4 



(16£ 2 + 28£ - 9) 3 ar csin ^ + ( 2 7 - 66£ + 320£ 2 - 240£ 3 + 64E 4 ) VI 



(A14) 



Finally, let u s thi nk of the DED of the generalized Plummer models corresponding to the DF given in equation J29l >. For this 
case, equation (IA5> reduces to a similar integral as equation d29i 



dM _ (27r) 5/2 r(l -/3)f E (E) 1 f 1 (1 - tf,P)(- 3 -V)/P- 
~dE~ 2F X F(3/2-/?) J E if/ 2 ®-® 



(iff - Ef' 1 -^ dt/f. 



(A15) 



Like the DF, equation ( IA15I reduces to a closed form function of E for particular values of p or /3. Among them are for/? =1/2 



dE J 2 



(2nf (1 - EPflv p + 1 



E 2 



(1 - E p y /P [(p + 2) - (2p + l)E p ] E p -\ 



(A16) 



fory3 = -1/2, 



dM =f(E) (27r)3 ( i - Ep y +4/p 

dE lE{ >4(p+4) 



£3 



1 4 
2 F 1 \l,- + l;-+2;l-E p 
P P 



D + 1 (l _ E p ) l/ P r n / 1 4 \ 

- — 7- t 1 — (P + 3)(p + 4) - (5/? 2 + 12p + 3)E p + 2p(2p + Y)E lp ^ _1 2 Fi 1, - + 1; - + 2; 1 - , (A17) 

p+44 L \ /? /? / 



and for p = 1, 
dM 



(2n) 5 ' 2 r(3 - 2®r(l -yg) (1 - £) 7/2 ^ 
d£ M } 2^-^(9/2-3^) £ 5 /2-/j 2 i 



2r(5 - 2yS)F(3 - 2yS) 
r(9/2-3yS)r(7/2-/3)' 



1 3 9 



1 



(1 - £) 7/2 - 3 Vi 1 - 2y6, 5 - 2# - -ft E \ 2 F l - - # - -ft - - 3ft 1 - E 



(A18) 



This suggests that dM/dE diverges as E — > for p < 1, that it vanishes for p > 1, and that it has a finite limiting value of 
8(2 —fi)/(5 - 2/3) for p = 1. We believe that t his is because the be havior of dM/dE near £ = is dominated by the stars at large 
radii and the r~ 4 fall-off marks a critical point (E vans & Anl2006l) . 



